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ABSTRACT 

We derive the radio spectral index of the Barnard Loop (BL) from large-scale radio 
surveys at four frequencies and find it to be a thermal source. We use the radio data 
together with Ha data to determine the electron temperature in BL, the A Ori HII 
region, and a high-latitude filament; all of these regions are somewhat cooler than 
typical HII regions. 

We perform least squares fits of the DIRBE diffuse IR intensities to the 21-cm line 
and radio continuum intensities. After the resolution of a "geometrical conundrum", 
this allows us to derive the electron density Ue] we find Ue ~ 2.0 cm~'^ and pressure 
f 24000 cm-3 K. 

Grains within BL are warmer than in HI regions. Trapped La accounts for the 
extra heating that is required. This is a general effect that needs to be accounted for 
in all analyses that examine IR emission from H"*" regions. Very small grains that emit 
60 radiation are enhanced in BL relative to HI by a factor of 2-3, while PAH's that 
emit 12 /im are probably deficient by a factor ~ 2. 



Subject headings: ISM: dust — HII regions — ISM: individual (Barnard Loop, A Ori) 



1. INTRODUCTION 

The Wisconsin Ha Mapper (WHAM) Sky Survey maps the Ha emission from the diffuse 
Warm Ionized Medium with unprecedented sensitivity and coverage (Haffner, Reynolds, and 
Tufte 1999). One of the brightest large regions is the Orion-Eridanus superbubble. We have been 
combining data from many sources for this region to gain a detailed understanding of physical 
conditions and processes (Heiles, Haffner, and Reynolds 1999). Here we present our results and 
analysis of a small portion of this region, the Barnard Loop (BL). 
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O'Dell, York, and Henize (1967) developed a coherent physical model based on UV, optical, 
and radio data. Their UV images are particularly interesting because the UV is dominated by 
scattered light and allows derivation of dust grain properties; they are also of historical interest, 
having been taken by astronauts aboard the Gemini 11 spacecraft. They model BL as an ellipsoid 
in which the volume density increases roughly as distance^, and suggest that the ionized portion 
of BL is surrounded by a neutral HI shell. They discuss the density structure as being produced 
by radiation pressure of light from the central stars pushing on the grains. Over a period 3 x 10^ 
yr this would produce the observed structure, which should have expansion velocity ~ 9 km s~^. 
We will defer discussion of the kinematics and correspondence with this model to a later paper. 

Here we restrict ourselves to a detailed analysis of radio, optical, and diffuse IR emission 
data to determine the physical conditions and dust grain properties in BL. Section ^ analyzes 
radio continuum surveys at four different frequencies and finds BL to be a thermal (optically 
thin free-free) emitter. Section ^ combines the radio and WHAM Ha data to derive upper limits 
on electron temperatures. Section ^ considers a discrepancy between the radio/Ha temperature 
and previous temperatures derived from the Nil and Ha lines and derives a correction factor 
for the absolute intensity scales of the radio and WHAM surveys. Section |^ adopts reddenings 
and turns the upper limits into actual temperatures. Section |^ combines the DIRBE IR, 21-cm 
line (Hartmann and Burton 1997), and radio continuum data to derive the BL electron density, 
grain temperature, and relative population of very small grains and PAH's. Finally, Section ^ 
summarizes our findings. 



2. RADIO EMISSION: THERMAL OR NONTHERMAL? 

Reich (1978) provides a very useful history of radio observations of the Barnard Loop (BL), 
which we will not repeat here. He concludes that there are significant background emission 
fiuctuations, and these prevent an accurate determination of the spectral index of the Loop itself. 
However, Haslam and Salter (1971), studied the bright top portion at 85.5, 240, and 408 MHz and 
found Tb cc zy-2-05±o.30^ consistent with a thermal spectrum. 

In hopes of narrowing this conclusion, we examined four surveys of radio continuum emission 
to determine whether BL is a thermal or nonthermal radio emitter. These were the 408 MHz 
survey of Haslam et al (1983), the 820 MHz survey of Berkhuijsen (1972), the 1420 MHz survey 
of Reich and Reich (1986), and the 2326 MHz survey of Jonas, Baart, and Nicolson (1998). The 
three lower-frequency datasets are available from http://www. mpifr-bonn. mpg. de/survey.hirn\ . 



Jonas provided us with the full-resolution version of the 2326 MHz survey. 

We were concerned that the surveys might exhibit scale errors. We addressed this problem by 
using the AOri HII region as a calibrator, assuming that it was an optically thin free-free (thermal) 
source with Tg cx v^"^-^. Using this as a calibrator has both good and bad points. For the good, it 
is large compared to any of the surveys' angular resolution and it is not so bright as to introduce 
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nonlinearities. However, its low intensity makes its derived flux susceptible to zerolevel and 
background intensity fluctuations, which are definitely present, particularly in the lower frequency 
surveys where nonthermal background radiation is brighter. For each survey we calculated the 
excess brightness of the HII region for four different choices of ON and OFF regions; our final 
choice used only the brightest part of the HII region for the ON and two immediately adjacent 
areas for the OFFs, as shown by the three areas outlined by dashed lines in Figure |. We adopted 
the 2326 MHz survey as a standard and derived the factors required to correct the other surveys 
to produce the thermal scaling. For the four choices, these factors ranged from 0.99 1.98, 
0.98 1.41, and 1.34 1.57 for the 408, 820, and 1420 MHz surveys, respectively; for our final 
choice illustrated in Figure |l|, the factors are 1.61 ± 0.50, 1.19 ± 0.22, and 1.45 ± 0.08. The quoted 
uncertainties are simply half the difference in ranges quoted above and should be overestimates, 
because the differences occur primarily because of background fiuctuations and our final choice 
has the OFFs closest to the ON. 

If the intensity scales of the surveys were identical, these factors would be unity. The first two 
have large uncertainties, but nevertheless seem to differ significantly from unity. The 1420 MHz 
survey certainly differs significantly from unity. Reich and Reich (1988) find that the 1420 MHz 
survey needs to be multiplied by 1.55 to agree with the 408 MHz survey; our derived value of 1.45 
differs from this factor by only 6%. 

We used these factors to make difference images. First we multiply each of the three low- 
frequency surveys by its factor. Then we multiply each survey image by the factor ( 2326'^MHz ) ' 
which converts all the temperatures to the 2326 MHz temperature if the emission is thermal. 
Finally, we convolve the 2326 MHz image to the same angular resolution and subtract it. The 
result is a difference image in which departures from the arbitrary and unknown zero-level offset 
represent nonthermal emission. If there is any morphological structure in the difference images 
that looks like the BL, then the Loop exhibits nonthermal emission. 

Figure |l| exhibits the three scaled difference images, together with the 2326 MHz image. The 
diagonal cutoffs in the images arise from the declination limits of the surveys. The dotted half-ring 
guides the eye to BL and the calibration rectangles illustrate offsets ±0.05 K from each difference 
image. All four images exhibit stripes; they run along the scanning direction and represent 
zerolevel errors from one scan to another in the original data. The three lower frequency difference 
images reveal artifacts and also background intensity fiuctuations. The 1420 MHz difference 
image, lower left, exhibits broad vertical stripes; these are zero level errors because they extend all 
the way to North celestial pole. The 408 MHz image exhibits fiuctuations over a broad angular 
scale, about | the size of the image, or about 8°; these are not as distinct at 820 MHz and might 
be partly or wholly spurious. 

The AOri HII region is invisible in the difference images of Figure |l|, which serves as a check 
on our method. BL is not recognizable in either the 820 or the 408 MHz difference image. The 
408 MHz image does exhibit excess emission ~ 0.12 K to upper right of center, but this feature 
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has no morphological relationship to BL. Moreover, there is no obvious evidence of this excess in 
the 820 MHz map, so the structure might be spurious. BL seems to be dark in the 1420 MHz 
difference image. However, the most noticeable dark portions lie on the dark broad stripes, so this 
appearance may be spurious. 

The calibration rectangles are easily recognizable at 820 and 1420 MHz, while BL is not. 
Thus, the nonthermal low-frequency emission from BL < 0.05 K. At 2326 MHz, BL's bright 
horizontal top portion (BL-TOP in Figure has Tg ~ 0.35 K and the left vertical portion 
(BL-LEFT) has Tb ~ 0.15 K at 2326 MHz; the hmit of ±0.05 K means that < 15% and < 33% of 
the emission is nonthermal at the lower frequencies. At 1420 MHz BL might exhibit a deficiency 
at this level, but this is contrary to nonthermal emission, which has a steeper spectrum than 
thermal emission. Therefore, these visual estimates of upper limits on nonthermal emission are 
quite conservative. 

One usually derives the spectral index of a morphologically distinct extended feature by 
making plots of intensities at many positions at one frequency versus those at a different frequency 
and performing a least squares fit. Instead, for each survey we fit the difference image dif f to 
the 2326 MHz brightness T2Z2&- For BL-TOP we obtained = 0.58 ± 0.25, 0.05 ± 0.01, 

and -0.25 ± 0.03 for the 408, 820, and 1420 MHz data, respectively; for BL-LEFT we obtained 
1.79 ± 0.03 and —0.34 ± 0.03 for the 408 and 1420 data, respectively. The significant negative 
slopes for the 1420 data probably reflect the coincidence of BL with the broad stripes, and the 
significant positive slope for BL-TOP at 408 is not revealed in the difference image except as a 
morphologically inconsistent bright blob. In our opinion, these least squares fits are meaningless 
because of artifacts and nonthermal emission fluctuations in the low frequency images. 

We conclude that BL exhibits no recognizable nonthermal emission. For BL-LEFT our 
conservative upper limit is that < 33% of the 408 MHz emission is nonthermal, and for BL-TOP 
< 15%. Nonthermal brightness temperatures typically oc z/"^''', so at 2326 MHz these limits 
become < 12% and < 4%, respectively. 



3. TR/Ha- THE ELECTRON TEMPERATURE FROM RADIO AND Ha 

EMISSION 

In this section we derive the electron temperature by comparing the 2326 MHz and Ha 
intensities, assuming that all of the 2326 MHz emission is thermal. The important quantity for 
these emission processes is the emission measure EM, 



EM = nil , (1) 

where rig is the total electron density. Electrons come from both H"*" and He+, so n^, = fu^+n{H^), 
where fHe+ accounts for the ionized He; below, we argue that fHe+ = 1- 
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Figure 1 compares the radio and optical images. To facilitate quantitative comparison of 
the radio and Ha emission we convert the usual intensity units of brightness temperature and 
Rayleighs to identical units that are close to the emission measure, normalizing to T4 = 1 and 
E{B -V) = 0: 

EM2326 = 1.80 X 10^TB,obs , (2a) 



EM Ha = 2.7AlHa,obs , (2b) 

where EM is in cm~^ pc, T4 the electron temperature in units of 10^ K, and E{B — V) the 
reddening in magnitudes. These equations derive from Mezger and Henderson (1967) and Haffner 
et al (1999), respectively. Here Ts^ohs is the observed 2326 MHz brightness temperature and 

a,obs the observed integrated line intensity in Rayleighs. These differ from the true intensities 
for two reasons. One is the error in absolute calibration; we define the intensity scale errors by 

TB,true = /2326?B,obs (3a) 
lHoi,tTue = fHalHa,obs (3b) 



The other is angular size, which is discussed below. With these definitions, we have 

EM = T0-35/2326^M2326 , (4a) 
EM = TrfHe-^e'-''''^^-''^fHaEMHa , (4b) 

Here we have used the extinction curve of Savage and Mathis (1979). Internal extinction, 
which occurs within the emitting region itself, is less effective than external extinction because 
a significant fraction of the extinction is scattering, and internal extinction does not scatter 
photons away from the observer (Mathis 1983). We have assumed that the extinction occurs in 
a foreground cloud sufficiently far from the emitting region that scattered Ha photons are lost, 
so that the standard extinction law applies. If our assumption is invalid, then reddening has less 
effect than above. Combining equations § we have 

^ _ f f2326EM2326 \ ^'^^ ^-4.25E(B-V) fr:\ 

J-A,R/Ha - -f 7 FIT \ e ^ ' . [b) 

KlHe+JHaEMHa/ 

Here the subscript R/ Ha means that the temperature is determined by combining radio and Ha 
data. 
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Figure |2| exhibits two sets of the 2326 MHz radio and Ha optical images side-by-side; one set 
exhibits BL and the other set a weak filament. We delineate four regions with dashed lines: the 
AOri HII region; BL-TOP, the top bright portion of BL; BL-LEFT, the left dim portion of BL; 
and FILAMENT, a weak filament lying near b = —40°. For each region we make a least squares 
fit for coefficients A and B in the equation 



^2326 = A + B- EMffa ■ (6) 

Table |l] presents the results and Figure ^ exhibits plots of EM2326 vs EMua for the four 
regions. Table |^ also presents the deduced temperatures T^^^n/Ha after multiplying the slopes 
B for the BL and FILAMENT entries by fs = 1.032 as described immediately below. It is 
not appropriate to multiply A by fs because A is an offset whose value is determined by the 
background, whose angular scale exceeds 7°. 

Equations |3| are not quite correct because of inadequate angular resolution, otherwise known 
as "beam dilution" . The optical data integrate over a l°-diameter circle and are spaced by about 
1°; BL and the AOri HII region are larger than 1° in all dimensions, so no correction is required. 
However, for objects < 7° diameter the 2326 MHz intensities must be corrected upwards by the 
factor 



fs = ^, (7) 

where r/ is the beam efficiency, Qs the source solid angle, and the solid angle of a 7°-diameter 
circle (Jonas et al 1998). The AOri HII region has diameter ^ 7° so there is no correction. Using 
the high-resolution Ha image of Isobe (1973), we model BL as a rectangle of dimensions 1.3° x 7°, 

and we take fs = ( J'"^ ) = 1.032, where ^1.3 is the solid angle of a 1.3° diameter circle. 



4. Tji/Ha- CORRECTION FOR ERRORS IN ABSOLUTE CALIBRATION 

The accuracy of Tjiju^ depends on the accuracy of the absolute calibrations of the 2326 
MHz and Ha datasets. Absolute calibrations are notoriously difficult so we cannot rely on them. 
Rather, we use other temperature determinations to adjust our results. 

Here we consider other temperatures determined by two techniques, both of which compare 
the [Nil] A6583 and Ha lines. One uses the line intensity ratio and one the difference in line width; 
we denote these temperatures by T/v/z/Ha and TAr/z^id/i^a^irf, respectively. Comparing these lines 
has the great virtue that they are measured with the same instrument, so systematic errors cancel. 
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4.1. A Two- Temperature Toy Model 

The only significant problem from comparing the [Nil] and Ha lines is possible temperature 
variations along the line of sight. The [Nil] line emissivity increases exponentially with 
temperature, while the Ha and radio emissivities decrease with temperature; thus, the [Nil] line 
preferentially samples high temperature regions. In a low-density HII region, the equilibrium 
temperature is independent of density. However, the temperature depends on distance from the 
star, becoming highest near the edge where the spectrum of ionizing radiation hardens. For 
example, models 3 and 4 by Rubin (1968) are not too dissimilar from the A Ori situation and 
exhibit variations T4 ~ 0.45 — > 0.80. We can get a rough idea of these effects by considering a 
two-temperature toy model, r4 = 0.45 and 0.80, with equal emission measures at each temperature. 
We show the results in Table ^. 

For all entries in Table ^, T^^nii/Ho = 0.630 and T^ j^/jj^ = 0.570. As expected, 
Tnii/Ho > Tji/jjai but by a modest amount compared to the difference between the physical 
temperatures of the two regions. In contrast, Tj^ji^ifijuawid depends very sensitively on the 
relative turbulent velocities in the two regions and can easily lie outside the physical temperature 
extremes that exist. We conclude that Tj^iiwid/Hawid must be used with caution, unless one can 
be sure that the nonthermal velocity is constant in the region. 

4.2. Temperature Comparisons for Two Regions 

Both [Nil] and Ha lines have been measured for two of our objects, the A Ori HII region and 

BL. 



4.2.1. XOri HII Region 

This HII region is excited by the star the 08III star A Ori; such a star has effective 
temperature 34000 K (Binney and Merrifield (1998). The radius of the He"*" ionization zone is 
about 0.35 the H"*" zone (Osterbrock 1989), so the He"*" volume is only about 0.04 the H+ volume; 
accordingly, we neglect the He"*" and take fHe+ = 1- 

Our upper limit is T^^R/Ha < 0-71 ± 0.05. The exciting star has reddening E{B — V) = 0.12 
mag (Diplas and Savage 1994). If all of this extinction occurs in front of the HII region, then we 
have T^^R/Ha = 0.43 ± 0.03. 

Reynolds and Ogden (1982) measured the [Nil] and Ha lines; applying equation (11) of 
Haffner, Reynolds, and Tufte (1999) gives T^^Nii/Ha = 0-57 ± 0.03. Also, Reynolds and Ogden 
found T^^Niiwid/Hawid = 0-60 ± 0.17 from a linewidth comparison; if the uncertainty were smaller, 
then the agreement would suggest that the temperature doesn't change much along the line of 
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sight. We adopt the temperature from the intensity ratio, T^^nii/Ho = 0-57 ± 0.03. 

Ti^Nii/Ha ~ Ti^R/Ha = 0.14 zb 0.04. The difference is far larger than the uncertainty. However, 
the results for the toy model in Table |2| suggest that we should expect a difference of a few 
hundredths]]. Suppose that the true value of T^^r/ho should be a few hundredths smaller than 
T4,Nii/Ha, say 0.54. 

There are two ways to make our derived T4^^R/Ha equal to 0.54. For one way, the effective 
E{B — V) is not 0.12 mag, but instead the value required to reduce our upper limit of 0.71 to 0.54; 
this would make the effective E(B — V) = 0.06 mag. This could occur if most of the extinction 
to the star A Ori occurs within the HII region itself. This is highly unlikely: the A Ori HII region 
has a total column density to its center N^, ^ 2 x W^^ cm~^ (Reich 1978), which corresponds to 
E{B — ~ 0.03 mag for the normal gas/extinction ratio (Bohlin, Savage, and Drake 1978). This 
is much smaller than the total reddening to the star. 

For the other way, we adjust the intensity scale errors to attain agreement. From equation ^ 
we require (fc)' '' = Miif , or ^ = 1.13 ± 0.06. 



4.2.2. BL-TOP 

Again we take f}je+ — 1 because BL is far from the exciting stars. Our upper limit 
is T^ R/jja < 0-49 ± 0.07. There is a previous radio determination by Gaylard (1984), who 
measured radio recombination lines and continuum intensities at three positions; the average of 
his temperatures is T^ rjr = 0.52 it 0.08. The agreement between our upper limit and his value is 
good, which implies that the reddening to BL-TOP is very low. 

Reynolds and Ogden (1979) measured the [Nil] and Ha lines; again, applying equation (11) 
of Haffner, Reynolds, and Tufte (1999) gives T^^^nii/Ho = 0-61 ± 0.03. Also, Reynolds and Ogden 
found T^j^jj^i^/jjawid = 0-76 ± 0.09 from a linewidth comparison. These two temperatures differ 



by 2a. From our discussion in Section ^ above, we eliminate T4^^Niiwid/Hawid from consideration 
and adopt the temperature from the intensity ratio, T^^^Nii/Ha = 0-61 ± 0.03. 

'^i,Nii/Ha ~ '^A,R/Ha = 0-12 ± 0.08. The difference is significant at the 1.5cj level. However, 
it seems that there is temperature structure within BL, so the expected difference is perhaps 
~ 0.02. Thus, we ascribe the larger observed difference to a scale error in T4^^R/Ha ^-nd suppose 
that the true value of T^ r/j^^ should be ~ 0.59. If there is no extinction, then we must adjust the 

intensity scale errors to attain agreement. From equation |^ we require ("^^^ = o 49^0 07 ' o^' 
^ = 1.11 ±0.09. 

JHa 



^The difference for the toy model is 0.06, but the toy model considers only two extreme temperatures instead of 
a continuous distribution and therefore overestimates the difference. 
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4.3. The Final Correction Factor 

Above in Section |4.2| we found two independent values for the ratio of radio to optical 
correction factors. They agree well even though their uncertainties are rather large. We adopt the 
unweighted average J^s^e _ ^ -^2 ± 0.07; the uncertainty does not include the errors in our guesses 

JHa 

for the proper differences T^^Nii/Ha — T4,R/Ha- 

We provisionally assign all of the calibration error to fija, for several reasons. Firstly, the 
present Ha data are not corrected for sky transmission variations, either night-to-night or airmass. 
Secondly, the Ha absolute calibration is ultimately tied to the brightness of a 1° diameter region in 
NGC7000 (the North American Nebula), and the absolute intensity of this region is uncertain at 
the level ^ 12%. Finally, Jonas (1999) has determined the uncertainty in the 2326 MHz intensity 
to be /2326 = 1-00^. 

For the remainder of the present paper, we adopt /2326 = 1 and fna = 0.89. 

For the ensuing discussion we do not use T^^R/Ha- Rather, we use values corrected by the 

factor y^yff^J — 1-23 and denote these corrected values by the symbol T^.R/Halcorr- These 
values are listed in Table ||. 



5- Tji/H^: DISCUSSION FOR INDIVIDUAL OBJECTS 

Equation |5| shows that the values for T^,R/Ha\corr listed in Table |^ are upper limits to the 
electron temperature because the reddening E{B — V) cannot be negative. 



5.1. The AOri HII Region 



Following our discussion above in Section [4.2.1 , we assume that all of the extinction occurs 



in front of the HII region and adopt E(B — V) = 0.12. For this choice, our upper limit becomes 
T4,^R/Ha\corr = 0.53 ± 0.04. By design, this is consistent with T4^^Nii/Ha = 0.57 ± 0.03 (Section Q). 

This is considerably smaller than typical temperatures of HII regions near the Sun. In 
particular, the Ori A and Ori B HII regions have measurements of T^^r/r ranging from 
~ 0.71 — > 0.86 (Here the subscript R/R means determined from the ratio of radio recombination 
lines to radio continuum; see Reifenstein et al 1970, Shaver et al 1983). It is also smaller than 
the predicted temperatures: for "standard" abundances, Osterbrock's (1989) Figure 3.2 should 
be reasonably representative of the A Ori HII region and predicts an equilibrium temperature of 
T4 = 0.70. We now discuss this discrepancy between observed and predicted temperatures for the 
A Ori HII region. 



The temperature of an HII region is governed primarily by the abundances of Nitrogen and 
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Oxygen. Lowering the temperature from the predicted T4 = 0.70 to T4 = 0.57 requires that the 
coohng be increased by a factor of 1.3. To attain the factor 1.3, the [Oil] and/or [Nil] abundance 
must increase by the factor 1.8. It is very doubtful that this increase occurs because of an increased 
gas-phase abundance produced by grain destruction, because Meyer et al (1997, 1998) find that 
only ~ 30% of the O and ~ 0% of the N are locked up in grains. 

An alternative to the observed temperature being as low as 5700 K is the fraction [N'^/N] 
being lower than unity, because the observed value T^^^Nii/Ha = 0.57±0.03 relies on the assumption 
that N+/N = 11+ /H (Haffner et al 1999). In fact, models by Sembach et al (1999) suggest that 
this assumption is incorrect; instead, N+/N = 0.7 II+/H. If this is correct, then the observed 
T^^Nii/Ha becomes 0.62 ± 0.03. This raises the temperature of the A Ori HII region and makes it 
closer to, but still different from, the theoretical prediction. 

We do not understand the causes for the smaller temperature of the A Ori HII region. 



5.2. BL-TOP 



Following our discussion above in Section 4.2.2 , we adopt E{B — V) = 0. This makes the 
temperature equal to our upper limit, so T^ j^/jj^lcorr = 0.61 =b 0.09. Again, by design this is close 



to T, 



i,NII/Ha 



0.61 ±0.03. 



5.3. BL-LEFT 

Our upper limit is T^ j^/jj^lcorr < 1-06 =b 0.18. This is considerably higher than BL-TOP. We 
expect the physical conditions in BL-LEFT and BL-TOP to be similar because they are part of 
the same structure; this implies that the true temperature is lower than our raw measurements 
imply and that extinction is important. If BL-LEFT has the same temperature as BL-TOP, then 
equation ^ predicts E{B — V) = 0.13 mag. 

This set of conditions is entirely in line with the ratio of Ha to H/3 line intensities, which we 
denote by Ra/3- Isobe (1978) presents a map of this ratio from which we estimate 



Ral3, BL-TOP 



0.80 , (8) 



Ral3,BL-LEFT 

and a good approximation to the ratio of line intensities (Martin 1988) is 

Raf3,theory ~ 2.^Kr^^'^ e^'^^^'^^~^^ . (9) 

These, together with the Table H's ratio of the slopes B for the two regions. 
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provide an independent estimate of the quantities Tratio = t4 bl left ^^'^ 

Ediff = [E{B - V)bl-top - E{B - V)bl-left]; we obtain Tratio ~ 0.58 and Ediff ~ -0.19 

mag, which agrees quite weh with our above results. 

We conclude that BL-LEFT has the same temperature as BL-TOP, T^ ^ju^ ~ 0.61. Its 
reddening E{B — V)k, 0.13 mag corresponds to a foreground Nh ~ 8.0 x 10'^'^ cm~^. 



5.4. FILAMENT 

Our upper limit is T^^^Rjuo^corr < 2.07 ± 0.75. The uncertainty is large, but an eyeball 
examination of Figure |3| persuades us that this value is reliable and is statistically significantly 
higher than in the other regions. 

One spot in this region, {a, 6) = (04'*00'", 2.0°), was observed by Reynolds and Ogden (1979, 
entry 3 in Table 1); these data give T^^^Nii/Ha = 0.55 ± 0.03. Applying equation ^, we find 
E{B — y) ~ 0.31 mag. This corresponds to Nh ~ 1.9 x 10^^ cm~^, which in turn corresponds to a 
100 fim brightness of 28 MJy ster~^ for the usual IR to HI conversion of Schlegel, Finkbeiner, and 
Davis (1998). The observed 100 ^m brightnesses are smaller than this, varying from ~ 6 ^ 13 
MJy ster~^. This seeming discrepancy is not meaningful because the HI column density is large 
and, moreover, the IR emission per H-nucleus ~ 0.7 — > 1.1 MJy ster~^ cmgg^, with the higher 
ratios correlated with the higher 100 fim brightnesses, which implies the presence of H2. Under 
these dense conditions, the true column density exceeds that indicated by the IR emission as 
discussed in section ^. 

We adopt the reddening required to make T^ ji/jja match T^^nii, 0.31 mag. This weakens the 
observed brightness by a factor --^3.7. From Figure ^, we see that the observed Ha brightness is 
EMh a ~ 90 cm ^ pc; thus, the unabsorbed brightness ~ EMjjq^ ~ 340 cm ^ pc. This makes the 
intrinsic optical brightness of FILAMENT 0.6 the optical brightness of BL. 



6. DIRBE FITS 

IR emission traces the total warm/cold column density [N{Htot) = N{HI)+2N{H2)+N{H+)] 
and 21-cm line emission traces the HI column density N(HI); thus the appropriately-scaled 
difference traces [2N{H2) + N{H~^)]. For H^, note the important difference between its IR and 
Ha (or radio) emission: IR traces column density N{H^), while Ha traces emission measure 
EM = N{H^)ne. This means that the comparison provides rig. 

The important scaling factor is B, the IR emission per hydrogen nucleus. The IR emission 
peaks near 100 fi. Global fits at 100 /im typically obtain B in the vicinity of 0.65 MJy ster~^ '^''^20 
(Schlegel et al 1998; Reach, WaU, and Odegard 1998; Heiles, Haffner, and Reynolds 1999). Here 
we adopt the somewhat lower value 0.62 of Arendt et al (1998) because they obtained the full 
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DIRBE spectral coverage. Let the emission per nucleus in ionized gas be bjj+B. Then, neglecting 
H2, we have 



IRobs = A + BN{HI) + bH+BN{H+). (10) 

To treat observational data, consider a least-squares fit for the the coefficients A, B, and C in 
the equation 



IRobs = A + BN{HI) + CneN{H+) . (11) 

Here IRobs is the DIRBE data, N{HI) is the integrated 21-cm line intensity, and N{H~^)ne is 
from either the 2326 MHz data or the velocity-integrated Ha line intensity. Comparing the last 
term in these equations gives 



bH+B 

ne = ■ (1. 

The radio data are better than the optical data for these fits because they have higher angular 
resolution and are unaffected by extinction; their higher noise is unimportant because BL is 
intense. The radio and optical data give comparable results. In actual practice we do not fit 
equation ^ but instead 



IRobs = A + BN{HI) + C'EM2326 • (13) 



From equation 4a it is clear that 

ne = T^-'''fHl+f2326bH+§ = OMbn^^ . (14) 

where the numerical value is for T^^R/Halcorr = 0.61, fHe+ = 1) and /2326 = 1- 

Reliable results require choosing appropriate regions, specifically ones with small N{HI) and 
no H2. High N{HI) regions are unsuitable for three reasons: (1) their possible saturation of the 
21-cm line makes it an invalid tracer of N{HI); (2) their associated extinction shields the interior 
from starlight, making dust grains cooler and reducing the IR emission per H atom, and (3) their 
probable associated H2 provides extra dust unrelated to HI, increasing the IR emission per H 
atom. 
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6.1. Electron Density and Grain Emission Spectrum 

Most of BL does not satisfy the low N{HI) criterion because there are very dense background 
molecular clouds. Nevertheless, we were able to locate two small portions of BL and its environs 
where molecular clouds are absent. These regions are outlined by dashed lines in Figure We 
performed fits at the six relevant DIRBE wavelengths. 

Figure ^ exhibits the logarithm of the IR emission in MJy ster~^ per N{HI)2o, where the 
subscript means units of 10^0 cm ^ ; this is equal to log B in equation |l^. This figure also shows 
the global average spectrum determined by Arendt et al (1998). Our spectrum agrees well with 
the global one. 



At 100 ^m, equation 14 gives = 0.84 it 0.15 and 1.16 it 0.09 cm ^ for the top and bottom 



regions, respectively. These differ by 1.8cj, which we regard as fair agreement, and we adopt the 
average jp^ = 1.0 cm~^. 

°H+ 

The ratio ^ is the ratio of the H"*" and HI IR spectra; if this ratio depends on wavelength, 
then the grain emission spectrum differs for ionized and neutral gas. Figure ^ exhibits this ratio as 
normalized to unity at A = 100 fj,m. The H"*" spectrum differs significantly from the HI spectrum. 
Relative to A = 100 //m, the 60 fim points are definitely larger, and the 12 and 240 //m points are 
probably smaller (the errors are large, particularly at 240 /im). 

Our least-squares fit allows us to predict IRpr = A + BN{HI) + C'E M2326- Thus the 
difference IRobs — IRpr (in least-squares terminology, the residual) is an approximate measure of 
N(H2) combined with the deviation of rig from its average, i.e. 



R — IRobs — IRpr — B 



2iV(i/.) + A'(iJ+)(l-^) 



(15) 



There is also a minor contribution from saturation of the 21-cm line. 



The lower right-hand quadrant in Figure ^ maps this difference, which depends on two 
terms — the "molecular" and "ionized" terms. The ionized term can be positive or negative: it is 
negative in regions of large Ue, where a small column density gives bright 2326 MHz emission. 
Thus, dense HII regions would stand out as deficiencies in Figure ^ if they were not associated 
with molecular clouds. However, the dominant contribution to this image is clearly the Orion 
molecular clouds, which produce the large, prominent burned-out bright blob. 



6.2. The Geometrical Conundrum 

Above we adopted = 1.0 cm^'^. The 2326 MHz brightness temperature ~ 0.13 K, which 
corresponds to N{H~^)ne = 160 cm^® pc. Combining these with bi{+ = 1 gives N{H~^)2o = 5.0 
and path length L = 160 pc. 
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This path length is very long: BL is at the distance of Orion, about 450 pc, so our 160 pc 
line-of-sight length is 36% of the total distance. This, in turn, implies that BL is a cylinder that 
happens to have its axis accurately pointed towards us. This is not only deeply unsatisfying, but 
also at odds with what we would expect from the simple geometrical and physical model of a 
thick-walled spherical shell. 

In the vicinity of the fitted regions, BL has radius R ~ 6.1° and apparent thickness T ~ 1.3°; 
these translate to lengths i? = 48 and T = 10 pc. Such a shell has a maximum tangential path 
length through its edge of Lgp^ ~ 60 pc. Our derived path length exceeds this by a factor of 2.7; 
not only that, it exceeds even the diameter by a factor of 1.6! This discrepancy is the geometrical 
conundrum. 

Clumping makes the conundrum even worse. Suppose that the line of sight within BL contains 
clumps that occupy a fraction / of the sightline. This does not affect any of our derived physical 
parameters, including L; however in this case L refers to the length of the line of sight within 
the clumps. The total length within which the clumps lie is L/f. This makes the geometrical 
conundrum worse by a factor of /. 

To resolve the conundrum comfortably, we need to decrease L to something less than its 
maximum of 60 pc. For discussion purposes, we will assume L = 40 pc. This means we need to 
decrease our derived path length by a factor of 4. 

Suppose that ionized regions exhibit more grain emission per N{Htot) than neutral regions; 
this means bH+ > 1. This increases rig by the factor bu+ and decreases L by the square of bH+. 
To resolve the conundrum in a comfortable way, we need 6^+ ~ 4 or bu+ ~ 2. 

We adopt this as the resolution to the conundrum. Thus, we adopt Ue = 2.0 cm^'^. BL has 
T4,R/Ha\corr ~ 0.61; thus, P4 = = 2neT4 K 2.4. 



THE BL BIG GRAINS ARE WARMER THAN HI BIG GRAINS 



Above in Section 3.2, we found the geometrical conundrum requires more 100 fim IR emission 



per grain than in HI regions, and we adopted b{j+ = 2. And in Section 6.1, we found the IR 
spectrum of BL to be different from that of HI. Both the increased 100 iim emission and the 
changed spectrum suggest that the BL grains are warmer than the HI grains. In the remainder of 
this section, we will explore this temperature increase. We shall assume that the number of 100 
/xm emitting grains per H- nucleus is identical in BL and HI regions. 

We perform least square fits of the IR spectrum to determine the big grain temperature 
Tbg- In all fits we use a z/^ emissivity law and weight the points by the inverse square of the 
uncertainties. In HI regions, the A ^ 100 /im radiation comes almost exclusively from big grains 
{BG) and 60 /xm radiation from very small grains {VSG; Desert, Boulanger, and Puget 1990). 
These are two distinct grain populations and in solving for the temperature of A = 100 /xm 
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emitting grains we must consider only the BG's, i.e. we must use IR intensities for A > 100 fxm. 
For the three HI grain spectra in Figure ^ we obtain Tbg = (18.2, 15.4, 17.5) K for the (global, 
top, and bottom) regions, respectively. These refer to HI regions. The uncertainties on our derived 
temperatures (top, bottom) are large, so below we will use the global Tbg = 18.2 K as a reference. 

In BL, Figure ^ shows 60 /im is much stronger than in HI regions by a factor of 1.9 to 3.4. 
This can occur for two reasons. One, the population of VSG's is enhanced in BL relative to HI 
regions. Two, the temperature of BG's is considerably larger in BL than in HI regions. 

Below, we will least-squares fit the grain temperatures to the IR spectra for the H"*" in BL 
(top, bottom) and exhibit the results in Table ^ The table also displays two important related 
quantities. One is bjj+, the ratio of the 100 fj,m emissivity at the derived Tbg in BL to that in the 

global HI. The other is ( Tbg'^hi ) ' ^ grain emissivity oc u"^, this is the expected ratio of grain 
heating (equal to total grain IR emission). 

To begin, we suppose that the enhanced 60 /xm intensity in BL comes from BG's that 
are warm enough to overwhelm the VSG emission. Therefore, in this fit we derive Tbg by 
including the amount of 60 fj.m intensity that exceeds what we expect from the the VSG's, 
i.e. the 60 //m emission above the dotted line in Figure ^. For H"*" in BL (top, bottom) we obtain 
Tbg = (24.2, 22.4) K. These correspond to bB+ = (6.9, 4.2), which is much larger than the = 2 
estimate above. This can be reconciled by making the clumping factor / 7^ 1. If the physical path 
length over which the H"*" is distributed is fixed at 40 pc, then / oc &]^+; if &_h-+ = 2, as we argue 

above, then we require / = (^■^^■'^^ = (11.9,4.4) and the path lengths actually occupied by 
the H+ in BL are L ~ (4, 9) pc. This is a high degree of clumping, and we would expect to see 
considerable small-scale structure of the Ha emission within BL. However, the image presented by 
Isobe (1973) doesn't give this impression. We conclude that the enhanced 60 /im intensity in BL 
does not come from warm BG's. 

We next consider the other alternative, namely that the enhanced 60 //m comes from an 
enhanced population of VSG's. In this fit, we derive Tbg in the usual way, by including only the 
A > 100 /im data. For H+ in BL (top, bottom) we obtain Tbg = (24.2,22.4) K. These correspond 
to bjj+ = (1.7,2.2), which agrees well with the bff+ = 2 estimate above. These grains require an 
excess heating rate relative to global HI by factors of (1.5, 1.9). 

This excess heating rate is very close to that expected from the added heating by La 
photons trapped in ionized gas. Spitzer (1978) shows that this mean intensity of these photons 
is Ilo ~ 1-1 X 10~^ne erg cm~^ s~^ str~^. The ISRF is equivalent to a blackbody at 3.14 K 
(Mathis, Mezger, and Panagia 1983), or Ijsrf ~ 1-8 x 10"'^ erg cm~^ s~^ str~^. For rig = 2 cm~^ 
(Section |6.2| ), we have ^'^i^jg^p^°' — 2-3, which is close to the excess heating obtained in the above 
paragraph. 

We emphasize that this extra grain heating is a generally-occuring process in all regions 
that trap La photons (the "on-the-spot" case). The H"*" grain heating rate ~ (1 -|- 0.6ne) times 
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larger than the HI grain heating rate. This translates into an increased IR emission per grain. 

We much prefer this second alternative because the it is consistent with the bjj+ = 2 estimate 



obtained from our geometrical argument in Section 6.2. Also, the increase in Tbg is modest and 
is attained with something close to the expected increase in grain heating rate. We conclude that 
the 60 /im excess in the of BL comes from an enhanced population of VSG's. 

In contrast, the relative 12 //m intensity from BL is smaller than that from the global HI. 
This result is strongly suggested from Figure ^ but, given the uncertainties, not absolutely certain. 
The 12 /xm emission comes from a third population of grains; many workers, including Desert et al 
(1990), believe they are PAH's. We conclude that, with high probability, the population of PAH's 
relative to BG's in the H"*" gas of BL is smaller than in global HI. 

Our conclusions, then, are that the VSG's are more abundant and the PAH's less abundant 
in BL than in global HI. This is similar to but firmer than the results of Desert et al (1990) on the 
H+ gas in the California nebula. Clearly, it would be desirable to confirm this trend by studying 
other HII regions. 



8. CONCLUSIONS 

In Section ||, we derived the radio spectral index of the Barnard Loop (BL) from large-scale 
radio surveys at four frequencies. We eliminated possible scale errors by assuming the HII region A 
Ori to be a thermal (optically thin bremstrahlung) radio source. We found BL to have a thermal 
spectrum, too. 

Having found the radio emission of BL to be thermal, we could combine the radio and 
Ha line data in Section ^ to derive upper limits on the electron temperature (Table |^) for four 
regions: the A Ori HII region, two regions in BL, and a high-latitude filament in the wall of 
the Orion-Eridanus superbubble. In Section |^ we discussed two of these regions in detail and 
compared our temperatures with those previously obtained from Nil and Ha line ratios. Our 
currently-derived temperatures were lower than the previous ones. We developed a toy model to 
explore the effect of temperature structure along the line of sight, and found that one expects 
the radio/Ha temperatures to be somewhat smaller than the NII/Ha temperatures. However, 
most of the discrepancy is a result of inaccurate absolute intensity calilbrations. Provisionally, we 
accepted the 2326 MHz radio survey calibration as accurate (Jonas et al 1998) and ascribed all 
the error to the WHAM survey data, whose final calibration has not yet been done; the current 
WHAM intensities need to be multiplied by the factor 0.89. 

We discussed the four regions in some detail and adopted reddenings, which enabled us to 
derive actual electron temperatures (Table ||). Temperatures in all four regions are somewhat 
smaller than in the Orion nebula and other HII regions near the Sun. 

In Section ^ we performed least-squares fits of the DIRBE diffuse IR intensities to the 21-cm 
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line and radio continuum intensities. Our derived IR spectrum spans A = 12 ^ 240 fim and agrees 
well with the spectrum from global HI derived by Arendt et al (1998). In this fit, the ratio of the 
derived coefficients of the 21-cm line and radio continuum allow one to determine the electron 
density Ue- However, our derived Ue is small, requiring a very long path length through the Ha 
emitting region, which is inconsistent with the BL morphology and leads to the "geometrical 



conundrum" of Section 6.2. 



The conundrum can be resolved if the H"^ grains emit more IR per H-nucleus than HI grains 
do. This, together with the modified IR spectrum from the H"'' gas, is explained in Section |^ by a 
higher grain temperature in the H"*" than in the HI {Tbg in Table ^: HI is the first entry, H"*" the 
last two entries). The increased grain temperature in the H"*" gas agrees well with that expected 
from the extra heating produced by trapped La photons. The H"*" grain heating rate is generally 
higher than the HI heating rate by a factor (1 + 0.6ne). 

In BL, Ue ^ 2.0 cm-3 and T 6100 K. The pressure is ^ ^ 24000 cm'^ K. In addition, the 
H''' grains in BL exhibit excess 60 //m emission and deficient 12 /im emission, indicating that very 
small grains (VSG's) are more abundant and PAH's less abundant in BL than in the global HI. 



We thank D. Finkbeiner for much pleasurable and instructive consultation and for supplying 
destriped versions of the IR and 408 MHz datasets in a convenient form; J. Jonas for a copy of his 
thesis, the 2326 MHz dataset, and instructive discussions; P. and W. Reich for instruction and 
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on the radio properties of BL. WHAM is supported by National Science Foundation grant 
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Fig. 1. — Three scaled difference images, witfi tliermal radiation and an arbitrary zerolevel 
subtracted away; and tlie 2326 MHz image, as labeled. Each image is about 26° in size and b 
increases upwards. The two black circles are the Ori A and Ori B HII regions. The dotted line 
guides the eye to BL; the dashed lines enclose ON and OFF regions for the AOph HII region. The 
rectangles aid the eye in interpreting the greyscale calibration (see text). 

Fig. 2. — Two pairs of 2326 MHz radio and Ha optical images side-by-side; each image is about 
26° in size and b increases upwards. The upper set exhibits BL and the AOri HII region; the two 
burned-out blobs are Ori A and Ori B. The bottom set exhibits the FILAMENT, which lies near 
b = —40°. The optical data integrate over a l°-diameter circle and are spaced by about 1°; the 
positions observed are marked by dots. The dashed lines enclose regions where least squares fits 
were done (see text). 

Fig. 3. — Data points and least squares fits for the four regions shown on Figure ^. 

Fig. 4. — Three images of data used for the DIRBE least-squares fits. The dashed lines outline the 
small regions where background molecular clouds are absent and allow a reasonably accurate least 
squares fit of equation |l^. Each image is about 26° in size and b increases upwards. 

Fig. 5. — B is the logarithm of the IR brightness (MJy ster~^) per 10^'^ H-nuclei in the neutral gas; 
the arrow marks the global average. Dashed and solid lines are for the top and bottom regions in 
Figure 1^, respectively; the dotted line and global average are from Arendt et al (1998). 

Fig. 6. — which is the ratio of the H"*" and HI IR spectra, normalized to unity at A = 100 ^m. 
Dashed and solid lines are for the top and bottom regions in Figure ^, respectively; the dotted line 
is the global average of Arendt et al (1998). If we had derived ng at wavelengths other than 100 
/im, the derived Ue would equal Mr times the 100 ^m value. 
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Table 1. ELECTRON TEMPERATURES 



REGION 


A 


B 


g-4.25E(B-V) 


T4,R/Ha\corr 

(upper limit) 


E{B - V) 
(adopted) 


T4^Ii/Ha\corr 

(adopted) 


AOri 


371 ±7 


0.83 ±0.03 


0.71 ±0.05 


< 0.88 ± 0.06 


0.12 


0.53 


BL-TOP 


350 ± 15 


0.67 ±0.05 


0.49 ±0.07 


< 0.61 ±0.09 





0.61 


BL-LEFT 


204 ± 17 


0.83 ±0.14 


0.86 ±0.14 


< 1.06 ±0.18 


0.13 


0.61 


FILAMENT 


161 ± 13 


1.63 ±0.59 


1.68 ±0.61 


< 2.07 ±0.75 


0.31 


0.55 



Note. — A and B are results of the least square fits of EM2326 to EMna (Equation q and Figure 
T4^ii/Ha\corr in column 5 is identical to column 4, but corrected for intensity scale errors (Section it is an 
upper limit. Column 6 is the adopted reddening. Column 7 is the actual value of T4^ji/jja\corr using the adopted 
reddening in Column 6. 



Table 2. TEMPERATURES DERIVED FROM A TWO-COMPONENT TOY MODEL 



VTURB, RGN 1 


VTURB, RGN 2 


'^4,NIIwid/ Hawid 


VTURB 


Hq width 


Nil WIDTH 


17 


17 


0.547 


17.19 


26.88 


21.39 


7 


7 


0.542 


7.47 


21.91 


14.75 


12 


7 


0.678 


7.86 


23.43 


15.10 


17 


7 


0.872 


7.96 


25.28 


15.36 


7 


12 


0.408 


11.72 


22.37 


17.17 


7 


17 


0.210 


15.85 


22.88 


20.05 



Note. — We assume regions 1 and 2 have equal emission measures with T4 — 0.45 and 0.80, respectively; and 
(f ) = 5.43 X 10"'\ For aU cases, = 0.229, Tr/hc = 0.570, and Tj^n/Ha = 0.630. Columns 3 and 4 are the 

temperature and nonthermal linewidth that would be derived from the observed Nil and Ha lines. All velocities 
are halfwidths in km s~^. 
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Table 3. GRAIN TEMPERATURES, 100 f^m EMISSIVITY, AND 

HEATING RATES 



CONDITION 


Tbg, K 


bH+ 




GLOBAL HI 


18.3 


1.0 


1.0 


H+ with 60 ^im, TOP 


24.2 


6.9 


5.4 


H+ with 60 nm, BOT 


22.4 


4.2 


3.4 


H+ w/o 60 //m, TOP 


19.6 


1.7 


1.5 


H+ w/o 60 ^m, BOT 


20.3 


2.2 


1.9 



Note. — GLOBAL HI temperatures are from our fit to the spectrum of Arendt et 

al (1998); all temperature fits assume a emissivity law, so Column 4 is proportional 
to the grain heating rate. In Column 3, 6^+ is the A = 100 jj,m grain emissivity in the 
H"*" relative to that in HI. 



EMISSION MEASURE EM 
128 255 382 510 




SURFACE BRIGHTNESS I 
0.0 6.2 12.5 18.8 25.0 




